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The heavy-Fermion liquid with short-range antiferromagnetic correlations is carefully considered in the two- 
dimensional Kondo-Heisenberg lattice model. As the ratio of the local Heisenberg superexchange Jh to the 
Kondo coupling Jk increases, Lifshitz transitions are anticipated, where the topology of the Fermi surface (FS) 
of the heavy quasiparticles changes from a hole-like circle to four kidney-like pockets centered around (tt, it). 
In-between these two limiting cases, a first-order quantum phase transition is identified at Jh / Jk = 0.1055 
where a small circle begins to emerge within the large deformed circle. When Jh /Jk — 0.1425, the two 
deformed circles intersect each other and then decompose into four kidney-like Fermi pockets via a second- 
order quantum phase transition. As Jh / Jk increases further, the Fermi pockets are shifted along the direction 
(7r, 7r) to (7r/2, 7t/2), and the resulting FS is consistent with the FS obtained recently using the quantum Monte 
Carlo cluster approach to the Kondo lattice system in the presence of the antiferrmagnetic order. 

PACS numbers: 64.70.Tg, 71.27.+a 



Quantum phase transitions are emergent phenomena ob- 
served in many strongly correlated electron systems and has 
attracted much interest. An electronic transition associated 
with the change of Fermi surface (FS) topology, the so-called 
Lifshitz transition[l], can be induced without any spontaneous 
symmetry breaking and local order parameter. The Lifshitz 
transition is assumed to be a quantum phase transition at 
T = and it becomes a crossover at finite temperatures. Elu- 
cidating the nature of the Lifshitz transition, it seems that the 
transition manifests itself dramatically[2, 3] only when other 
degrees of freedom like lattice or spin couple strongly with 
the electronic states. 

Heavy fermion materials have played a particularly impor- 
tant role in the study of quantum critical phenomena. The 
Kondo lattice model is believed to capture the basic physics 
of heavy fermions. The model describes a lattice of local 
spin- 1/2 magnetic moments coupled antiferromagnetically to 
a single band of conduction electrons. The huge mass en- 
hancement of the quasiparticles can be attributed to the co- 
herent superposition of individual Kondo screening clouds, 
and the resulting metallic state is characterized by a large FS 
with the Luttinger volume containing both conduction elec- 
trons and localized moments. Competing with the Kondo sin- 
glet formation, the localized spins indirectly interact with each 
other via magnetic polarization of the conduction electrons - 
the Ruderman-Kittel-Kasuya-Yosida interaction. Such inter- 
action dominates at low values of the Kondo exchange cou- 
pling and is the driving force for the antiferromagnetic ( AFM) 
long-range order quantum phase transitions [4, 5], In a recent 
experiment^] a jump in the Hall coefficient for YbRh2Si2 has 
been observed, and a sudden change in the FS topology from 
a large FS to a small one was suggested at the magnetic quan- 
tum critical point. The nature of this phase transition is cur- 
rently under hot debate[7-13]. 

So far most of investigations focus on the possible FS 
reconstruction around the magnetic quantum critical point. 
However, we would like to point out that the FS topology 



in the paramagnetic heavy-fermion liquid phase may also be 
drastically changed by the short-range AFM spin correlations 
between the localized spins, leading to the Lifshitz phase tran- 
sitions. The nature of such a quantum phase transition has not 
been thoroughly explored yet. It is well-known that the large- 
N fermionic approach can be used to treat the Kondo lattice 
model in the Kondo singlet regime very efficiently, leading to 
the paramagnetic heavy-Fermion liquid state[14-16]. To con- 
sider the effects of the short-range AFM spin correlations, it 
is more straightforward to explicitly introduce the Heisenberg 
AFM superexchange Jjj between the localized spins to the 
Kondo lattice system[8, 9, 12, 17, 18]. 

In this paper, we apply the large- N fermionic approach to 
the Kondo-Heisenberg model on a two-dimensional square 
lattice in the limit of Jk > Jh- By introducing uniform 
short-range AFM valence-bond and Kondo screening param- 
eters, a fermionic mean-field theory is carefully re-examined, 
and such a mean-field theory becomes exact when the degen- 
eracy of the localized spins N becomes infinite. Away from 
half-filling, at the conduction electron density n c = 0.9, for 
example, as x = Jh/Jk increases, we find that the topol- 
ogy of the FS of the heavy quasiparticles changes from one 
hole-like circle to four kidney-like pockets around (w, w). In- 
between these two distinct limits, we will identify a first-order 
quantum phase transition at x lc = 0.1055, where a small cir- 
cle begin to emerge within the large deformed circle. Then the 
inner circle gradually expands, deforming to a rotated squared 
circle. Whence = 0.1425, the two deformed circles intersect 
each other and then decompose into four kidney-like Fermi 
pockets, resulting in a second-order quantum phase transition. 

The model Hamiltonian of the Kondo-Heisenberg lattice 
model is given by: 

H = ^ fkCkcrCkcr + Jr ^ S s ; • Sj + J H ^ Sj • Sj, (1) 
k,er i (ij) 

where ct creates a conduction electron on an extended or- 
bital with wave vector k and z-component of spin a =t, I- 
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The spin- 1/2 operators of the local magnetic moments have 
the fermionic representation Sj = \ J2aa> fia T >?>?'fi<?' w i m a 
local constraint flcrfic = 1' where r is the Pauli matrices. 
Following the large- N fermionic approaches, 17] the Kondo 
spin exchange and Heisenberg superexchange terms can be 
expressed up to a chemical potential shift as 



•- V f f f- ^ f- , 

2 / j J i<j J J a J jcr'Jta ; 



■3<r c i<r- 
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then a uniform short-range AFM valence bond and Kondo 
screening order parameters can be introduced as 
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FIG. 1: (Color online) Schematic plot of the conduction electron 
band ek hybridizing with the very narrow band Xk in the n c = 1 
case. The dashed lines denote the original bands, while the solid 
lines correspond to the hybridized bands, (a) The short-range AFM 
order parameter x > 0, (b) \ < 0. Clear distinction of these two 
different cases shows that there is always a finite gap in case (a), 
while being gapless in case (b). 



Generally speaking, apart from the uniform short-range AFM 
valence bond state, there are other possible competing states, 
including various flux phases or plaquette states. However, 
the stabilization of those states would require the presence of 
the translational lattice symmetry breaking and/or additional 
frustrating interactions. Since there is no evidence of transla- 
tional symmetry breaking in the paramagnetic heavy-fermion 
liquid states, we will not consider those effects at this stage. 
Although most heavy fermion systems are three dimensional, 
as far as the FS topology is concerned, it is conceptually sim- 
pler in the model discussion to start with a two-dimensional 
case of this Kondo-Heisenberg lattice model. 

To avoid the incidental degeneracy of the conduction 
electron band on a square lattice, we choose ej. = 
—It (cos k x + cos k y ) + At' cos k x cos k y — /i, where t and t! 
are the first and second nearest neighbor hoping matrix ele- 
ments, respectively, while \i is the chemical potential, which 
should be determined self-consistently by the density of the 
conduction electrons n c . Under the uniform mean-field ap- 
proximation, the f-fermions/spinons form a very narrow band 
with the dispersion Xk = Jh X ( cos k x + cos k y ) + A, where 
A is the Lagrangian multiplier to be used to impose the local 
constraint on average. 

In general, there are two interesting mean-field phases. One 
is the uniform short-range AFM ordered phase in the limit 
of Jh > Jk where V — but ^ / 0. In this phase, the 
spinons represented by f-fermions are unconfined and have a 
dispersion Xk- Another phase is the heavy electron phase with 
short-range AFM spin correlations in the limit of Jh < Jk, 
where both V ^ and \ ^ 0. Actually the former limit has 
been extensively studied in many previous investigations,^, 
13, 17] but we will focus on the latter limit. 

Thus the corresponding mean-field Hamiltonian reads 
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with E = N (-A + Jhx 2 + JrV 2 /2). The quasiparticle 
excitation spectra can be easily obtained 

1 



_(±) 



[(e k + Xk) ± Wk] , 



(3) 



which implies that the conduction electron band e k has a fi- 
nite hybridization with the spinon band Xk- Here Wk = 



(e k — Xk) + (JrV) . Accordingly, the ground-state en- 
ergy density can be evaluated as 
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A + Jh X 



\JkV 2 , (4) 



where 8 (— s) is the theta function. Then the self-consistent 
equations for the mean-field variables \, V, an d A can be de- 
rived by minimizing the ground state energy = 0, ^7- = 







t)s u 



-g^- — 0, and the chemical potential \i should be deduced 
from the relation n c = — ■ This leads to the following 



self-consistent equations at zero temperature, 
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(5) 



where 7k = cos k x + cos k y . In the following, we will assume 
that t' /t = 0.3 and n c = 0.9, which is away from half-filling 
and in the paramagnetic metallic phase. [1 1] 

The most important thing is to notice that the resulting two 
renormalized quasiparticle bands crucially depend on the 
sign of the AFM order parameter x- If X is positive, e^ - ' 
and e^ + ' are separated by an indirect energy gap as displayed 
in Fig. la, where the conduction electrons hybridize with the 
hole-like f-fermions/spinons. For the negative value of x> 
however, the renormalized quasiparticle bands and 
always have a finite overlap as shown in Fig. lb, where the 
conduction electrons actually hybridize with the particle-like 
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FIG. 2: (Color online) The lower renormalized quasiparticle band 
in the direction (n/2,ir/2) —¥ (it, it) (n, 7r/2) as increasing 
the strength of the AFM spin fluctuations, (a) to (i) correspond to 
J K /t = 2.0 and J h /Jk = 0.05, 0.1055, 0.1056, 0.11, 0.125, 
0.1375, 0.1425, 0.145, and 0.25, respectively. In the calculation 
t' /t = 0.3 and n c = 0.9. The dash line denotes the Fermi level. 



FIG. 3: (Color online) The Fermi surface corresponding to the lower 
renormalized quasiparticle band shown in Fig. 2. The ratio Jh / Jk is 
indicated in graphes. Here the center of the Fermi surface has been 
shifted from (0, 0) to (tt, tt) in order to clearly exhibit the topology 
of the Fermi surface. 



f-fermions/spinons. Even in the half-filling case n c = 1, the 
model system has the properties of a semi-metal with a FS 
consisting of one electron-like pocket around k = (0, 0) and 
one hole-like pocket around k = (it, tt). Then the resulting 
ground state will exhibit an instability towards the AFM spin- 
density wave or s±-wave paring superconductivity. [19] 

When the self-consistent calculations are carefully per- 
formed, we find that the mean-field AFM order parameter \ 
is always positive in the range < Jh/Jk < 0.5 so that 
the resulting state is a stable paramagnetic metal. In Fig. 2, 
as the strength of the AFM spin fluctuations grows up, we 
present the evolution of the band structure of the renormal- 
ized heavy quasiparticles around the Fermi level in the di- 
rection (7r/2,7r/2) — > (7r,7r) — > (tt, 7r/2) of the first Bril- 
louin zone. Fig. 2a to Fig.2i correspond to Jx/t = 2.0 and 
Jh/Jk = 0.05, 0.1055, 0.1056, 0.11, 0.125, 0.1375, 0.1425, 
0.145, and 0.25, respectively. 

For a fixed JK/t, we can clearly see in Fig. 2a and Fig. 2b 
that the M point (tt, tt) is the maximum of the lower renormal- 
ized quasiparticle band e^™' for Jh/Jk < 0.1055, and then 
the Fermi level (ep = 0) crosses the quasiparticle band only 
once in the the direction (tt/2,tt/2) — > (tt,tt) or (tt,tt) — > 
(tt,tt/2). For Jh/Jk > 0.1055, the M point becomes a local 
minimum of the lower renormalized quasiparticle band ^ 
in Fig. 2c to Fig.2i. In the direction of (tt/2, 7t/2) — > (ir, tt), 
the Fermi level always crosses the quasiparticle band twice. 
However, in the the direction of (it, it) — > (tt,tt/2), the 
Fermi level crosses the quasiparticle band twice only when 
0.1055 < Jh/Jk < 0.1425, while the Fermi level can cross 
the quasiparticle band once for Jh/Jk > 0.1425. Therefore, 



both Jh/Jk = 0.1055 and Jh/Jk = 0.1425 represent two 
special coupling strengths. 

Once the renormalized quasiparticle band structure is avail- 
able, the corresponding FS can be easily obtained. Corre- 
sponding to the band structure shown in Fig. 2, the FS are 
displayed in Fig. 3, where we have shifted the center of the 
FS from (0,0) to (tt, tt). For a fixed Jx/t = 2.0, we can 
clearly see that the FS is a hole-like circle around (tt, tt) for 
the parameter range < Jn/t < 0.1, and then the shape 
of FS deforms to a square for 0.1 < Jh/Jk < 0.1055. At 
Jh/Jk = 0.1056, the topology of the FS starts to change: 
a small circle emerges in the center of the deformed large 
square FS. As Jh/Jk is further increased, both circles ex- 
pand and the small one is deformed into a rotated square. Up 
to Jh/Jk = 0.1425, the two deformed circles intersect each 
other and then decompose into four kidney-like Fermi pock- 
ets. When Jh/Jk continues to increase, the resulting FS 
(not included here) will be shifted outward along the direc- 
tion m — > r. 

Recently, a quantum Monte Carlo cluster approach has 
been proposed to study the evolution of the Fermi surface 
across the magnetic order-disorder transition in the two- 
dimensional Kondo lattice system[ll]. In the AFM long- 
range ordered phase, the Kondo screening does not break 
down, and the heavy fermion bands drop below the FS giv- 
ing way to hole pockets centered around k = (tt/2,tt/2) and 
equivalent points. These results are fully consistent with the 
FS obtained by our calculation in the range 0.25 < Jh/Jk < 
1. 

In order to study the topology changes of the FS from the 
quantum phase transition aspects, we calculate the ground 
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FIG. 4: (Color online) (a) The ground state energy density e g as a 
function of x — Jh/Jk- (b) The first-order derivative of e g with 
respect to the parameter x. (c) The ground-state phase diagram of 
the system. 
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FIG. 5: (Color online) The quasiparticle mass enhancement factor as 
a function of x — Jh / Jk ■ 



state energy density e g and its first-order derivative with re- 
spect to the ratio of the coupling parameters x = Jh / Jk ■ The 
numerical results are displayed in Fig.4a and Fig. 4b. We find 
that two non-analytical points appear in the ground state en- 
ergy density. e g is finite and continuous in the parameter range 
< x < 0.5, However, its first-order derivative has a large 
jump at x\c = 0.1055, corresponding to a first-order (discon- 
tinuous) quantum phase transition. Moreover, a small kink 
appears at x^ c = 0.1425 in the first-order derivative, which 
corresponds to a jump in the second-order derivative of e g . So 
X2c denotes a second-order (continuous) quantum phase tran- 
sition. Actually both quantum phase transitions belong to the 
category of Lifshitz phase transitions. The ground state phase 



diagram is delineated in Fig.4c, where there exist three differ- 
ent paramagnetic heavy-fermion liquid phases: the conven- 
tional heavy-fermion liquid in x < 0.105, the heavy-fermion 
liquid with strong AFM spin fluctuations in 0.1425 < x < 1, 
and the intermediate phase 0.105 < x < 0.1425. 

Moreover, the effective mass of the heavy quasiparticle ex- 
citations is a function of the band curvature, so the topological 
changes of the Fermi surface can be reflected in the effective 
mass, which is related to the non-interacting band mass by the 
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FS. Hence the mass enhancement factor is given by 
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which is displayed in Fig. 5. Here we can also observe the two 
successive quantum phase transitions at Jh/Jk = 0.1055 
and 0.1425, respectively, consistent with the results from the 
analysis of the ground state energy density. Actually, since 
the effective mass enhancement factor is related to the opti- 
cal conductivity in infrared spectroscopy measurements, the 
above Lifshitz phase transitions can be observed experimen- 
tally. 

In conclusion, we have carefully studied the heavy-fermion 
liquid state in the two-dimensional Kondo-Heisenberg lattice 
system. As Jh/Jk grows up, the topology of the quasipar- 
ticle FS rapidly changes from one hole-like circle in the con- 
ventional heavy-fermion liquid state to four kidney-like pock- 
ets centered around (ir, tt), which is very close to the FS near 
the AFM magnetic quantum critical point. Between these two 
distinct FSs, a first-order quantum phase transition occurs at 
Jh/Jk = 0.1055, where a small circle emerges within the 
large deformed circle. When Jh/Jk = 0.1425, the two de- 
formed circles intersect each other and then decompose into 
four kidney-like Fermi pockets, and a second-order quantum 
phase transition takes place. Both quantum phase transitions 
belong to the category Lifshitz phase transitions. 

To some extent our present mean-field theory captures the 
heavy-fermion liquid physics of the Kondo-Heisenberg lattice 
systems, especially the Fermi surface evolution of the renor- 
malized heavy quasiparticles as the short-range AFM spin cor- 
relations between the localized magnetic moments are gradu- 
ally increased. In order to put the present results on a more 
solid ground, further investigations including the gauge fluc- 
tuations associated with the mean-field order parameters are 
certainly needed. 
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